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Abstract. We study inflation with multiple vector fields. In the presence of non- 
trivial couplings between the inflaton and the vector fields, it turns out that no-hair 
conjecture does not hold and vector-hair appears. In the case of uniform couplings, 
nevertheless, we find that the universe approaches an isotropic final state after transient 
anisotropic inflationary phases. For general couplings, we numerically show attractors 
are anisotropic inflation. Even in these cases, it turns out that the inflation always 
tends to minimize the anisotropy in the expansion of the universe. 
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1. Introduction 

It is believed that inflation accounts for features of cosmic microwave background 
radiation (CMB) observed by WMAP pQ. The remarkable predictive power of inflation 
is associated with the cosmic no-hair conjecture which claims that any classical hair 
will disappear exponentially fast and the universe asymptotically approaches de Sitter 
spacetime once the vacuum energy dominates the universe. This conjecture is supported 
by the no-hair theorem which holds in the special situations [21 E]. Thus, remaining 
quantum vacuum fluctuations are responsible for seeds of the large scale structure of the 
universe. The nature of quantum fluctuations reflects symmetry of de Sitter spacetime 
and gives rise to the following robust predictions. Firstly, the power spectrum of 
quantum fluctuations should have the statistical isotropy due to the rotational symmetry 
of the spatial slices. Secondly, the isometries of de Sitter result in constant Hubble 
expansion rate and the existence of the preferred Banch-Davis vacuum, which in turn 
implies the scale independent power spectrum. Finally, the slow roll conditions are 
guaranteed by shift symmetry in field space which implies suppression of nonlinearity 
and hence Gaussian statistics of fluctuations. Indeed, these predictions have been 
observationally confirmed. 

It should be stressed that the cosmic no-hair conjecture is just a conjecture. 
Historically, there have been challenges to the cosmic no-hair conjecture [H El El [7J 
IH [91 [101 [TT]. Unfortunately, it turned out that these models suffer from either the 
instability (121 [131 E]? or a fi ne tuning problem, or a naturalness problem. Recently, 
however, stable anisotropic inflationary solutions in the context of supergravity have 
been found as counter examples to the cosmic no-hair conjecture [151 [16]. More precisely, 
it has been shown that, in the presence of a vector field coupled with the inflaton, there 
could be a small anisotropy in the expansion rate which never decays during inflation. 
Subsequently, the generality of anisotropic hair in inflationary universes has been also 
investigated [HI HH [191 1201 [2H [221 [231 [21] . 

In spite of the failure of the cosmic no-hair conjecture, since the anisotropic inflation 
is an attractor solution, the predictive power of the model still remains [251 EHJ [271 121] . 
The interesting point is that the small anisotropy results in observationally relevant 
signatures. In fact, there exist imprints of the anisotropic expansion in the CMB [29] . 
Therefore, it is worth studying to what extent the cosmic no-hair conjecture can be 
violated. To this end, we examine inflationary models with multiple vector fields. In 
the presence of non-trivial couplings between the inflaton and the vector fields, the 
energy density of vector fields never dies away, which should give rise to anisotropy in 
the cosmic expansion. However, it has been known that multiple vector fields can take a 
configuration for which the universe is nearly, or entirely isotropic [9]. So, the question 
is the fate of the cosmic anisotropy in the presence of non-vanishing vector fields. It 
turns out that the accelerating universes tend to minimize their spatial anisotropy, 
which might be regarded as a generalization of cosmic no-hair conjecture, i.e., the 
cosmic minimum-hair conjecture. For uniform gauge-kinetic coupling models, we find an 
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isotropically inflating spacetime is the attractor whenever there are three or more vector 
fields. Although the analysis is carried out for Abelian vector fields, it also suggests that 
non-Abelian gauge fields, which are described by multi-component vector fields, cannot 
be anisotropic in the asymptotic future if the inflation is eternal. In more general cases, 
even though the anisotropy will not disappear, it will be minimized during inflation. 

The organization of the paper is as follows. In section II, we define models and give 
basic equations. In section III, we show the existence of isotropic inflationary solutions 
and the stability of isotropic solutions. This highlights a tendency of isotropization of 
the inflating universe. In section IV, we examined the global structure of the phase 
space. In particular, we investigate the fate of the anisotropy for general couplings. The 
final section is devoted to a brief discussion on phenomenology and the conclusion. 



2. Models 



,5' 



In this section, we setup models. For simplicity, we take power-law type inflationary 
models for which we have exact solutions even in the presence of vector fields [17] . 
We consider the following action. 

J d 4 xV^h L^R - \d.&r* - Voe~ k * - \ f; e^F^F^A , (1) 

where R is the Ricci scalar calculated from the metric h ab and G is the Newton constant. 
The N copies of Abelian gauge field A a m ^ with field strength F^ = d a A^ — d b A a m ^ 
are coupled to an inflaton <p with coupling constants g m . We consider the exponential 
potential with parameters Vq and k. The Latin letters denote the full space-time indices. 
Greeks are reserved for spatial parts. 

Following Ellis and MacCallum [30], we introduce the unit normal u a of the 
homogeneous hypersurface in generic Bianchi I space-times. The evolution of the spatial 
slice is described by 

u a] b = H{h ab + u a u b ) + a ab , (2) 

where H is the averaged expansion rate and the trace-free tensor a ab measures anisotropy 
of the universe. The decomposition of energy-momentum tensor can be carried out as 

Tab = pu a u b + 2q (a u b) + p(h ab + u a u b ) + 7T ab . (3) 

p, p, q a an n ab are respectively energy density, pressure, energy flux and anisotropic 
pressure seen by the observer with 4-velocity u a . The electric fields seen by an observer 
with the four-velocity u a are defined by 

E M = F ( r V (4) 

and we assume, again for simplicity, that the magnetic parts vanish, that is 

VabcdF {m)hc u d = (5) 
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where r? a bcd is the volume four-form. Then the matter variables are given as follows: 

1 i N 

p =-<p 2 + V e- k * + e 9 ^El m) E {m)a (6) 

2 2 

m=l 

q a =o (7) 

1 i N 

p = -0 2 - V Q e~ k<t> + e^E^E^ (8) 

2 6 

m=l 

vr afe = e ^ (-£< ro) ^ m) + ^E^E^ c {h ab + . (9) 

Here overdots denote differentiation with respect to the proper time t associated with 
the unit normal u a . 

In a group invariant orthonormal frame, the Einstein equations take the form given 
in Wainwright and Ellis [21] • Let us assume that E [ a ] and E [ a ] are non- vanishing and 
not parallel to each other. We choose the spatial frame such that 

EV = (E, 0,0), £f = 0. (10) 

This fixes the gauge completely and results in 

fil = Cr 23 , Q, 2 — — <7l3, fi 3 =(Ti2, (11) 

where f2 a represents rotational velocity of the frame with respect to the Fermi 
propagated one. In order to obtain a dimensionless dynamical system, we carry out 
expansion normalization. We adopt the following conventions: 

E± = 2^(^22 ±<733), S ^ = ^f fOT a ^P> ( 12 ) 

and 

<f>-Z- yl) - V ° e F - C(2) _ e 2 ^" _ e 2 /ION 

s# 2 ' fc V6^' a " Vqh ' a Vqh ,[ } 

where capital Latins run from 3 to N. The convenient choice of time coordinate r is 

dr = Halt (14) 

and we use primes to denote the time derivatives. The dynamical system is defined by 
the followings. The evolution of geometry is governed by the following equations 

E' + = (q- 2)E + - (E 2 2 + E 2 3 ) + 28 2 + 2(£{ 2) ) 2 - (£ 2 (2) ) 2 

N 

+ ^(2(^) 2 -(4 A) ) 2 -(4 A) ) 2 ), (15) 

A=3 

N 



E'_ = (g - 2)E_ - E 2 2 + E 2 3 + 2E 2 3 - 3(£f) 2 - 3 £((5^) 2 ~ ^f) , (16) 

A=3 

TV 

E' 12 = (g - 2 + 3E+ + E_)E 12 + 2E 13 E23 - 6S[ 2) S^ 2) - 6 ^ £[ A) £ { 2 A) , (17) 



A=3 
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TV 

E' 13 = (q - 2 + 3S + - £_)£ 13 - 6j2 £ i A)£ z A) > ( 18 ) 

TV 

£23 = (? " 2 " 2S -) S 23 - 2S 12 S 13 - 6 E 2 A) 4 A) ■ (19) 

A=3 

The equations for scalar field read 

3 TV 

$' = (g - 2 )$ + 3M> + 3^£ 2 + 3^ 2 {(4 2) ) 2 + (ff) 2 } + 3<m £ $>< A) ) 2 , (20) 

a=l A=3 

¥' = (2g + 2-fc$)# . (21) 
The vector fields obey the following equations 

1 - !$ - 2£ + )£, (22) 





= (? 


(^ (2) : 


i' = (g 


(4 2) : 


i' = (<f 



and 



(23) 



L - |$ + S + + E_)£ 2 (2) , (24) 

= (q - 1 - ^ - 2^ + )S[ A) + 2£ 12 £ 2 (A) + 2£ 13 4 4) , (25) 
(4 A) )' = (q - 1 - + E + + E_)£f> + 2£ 23 4 4) , (26) 

(£ { 3 A) )' = (g - 1 - + s + - S_)4 A) . (27) 

There is a constraint among the variables 

1 = E 2 , + ^(S 2 + S 2 2 + S 2 3 + S 2 3 ) + h 2 + * 

3 TV 

+^ 2 + (^ 2) ) 2 + (4 2) ) 2 + EE(^ A) ) 2 - C 28 ) 

a=l A=3 

The deceleration parameter q = —1 — H/H 2 is given by 

9 1 3 N 

q = 2Z 2 + + Z -(^ + Z 2 12 + Z 2 13 + Z 2 23 ) + 

a=l A=3 

Now, we are in a position to study the fate of anisotropy in the inflation with 
multiple vector fields. 

3. Cases of uniform couplings g = g m 

In this section, we restrict ourselves to the cases of uniform coupling constants g = g m 
and show that isotropic inflation is an attractor in the phase space for the number 
of vector fields greater than two. This clearly illustrates the existence of a principle 
replacing the cosmic no-hair conjecture. 
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3.1. Multi-vector Isotropic Inflation 



The system of equations ( fl~6j) to ( l29l) forms a dynamical system of dimension 3(N + 1). 
To understand the solution space, it is essential to find out equilibrium points and their 

(2) 

linear stability. Let us look for equilibrium points which have non- vanishing E,E\ and 
E 



For their time derivatives to vanish, it is required that 
9, 



q 



o. 

,(2) 



Then, from the evolution equations for E^' and > we obtain 

S12 = S13 = S23 = 0. 

Let us define the overall energy density parameter of the electric fields 
E 2 = E 2 + (E[ 2) ) 2 + {E^f + Y,^ A) f- 



Using the equilibrium conditions for \l/ and $ yeilds 



k + g 



£< 



k(k-g) 



(30) 

(31) 
(32) 

(33) 



k-g' k-g (k-g) 2 

Thus the existance condition for this class of solutions is kik — g) > 4. In order to have 
an inflating (accelerating) universe, we also require q < k 2 — g 2 < 0. The remaining 
task is to determine the relative strengths and angles among the field vectors through 
= 0. There are six equations to satisfy for 3(iV — 1) variables. 



2 IE 2 + (E[y + J2& A) ) 2 = (4 2) ) 2 + £(^) 2 + £(*: 



(4 2) ) 2 + £(4 A) ) 2 = £(4 A) ) 2 , 

A A 

,(2) c (2) 



(A) C (A) 
2 



Y^E^E 



(A) MA) 



(A) P (A) 
3 



0. 
0. 



E 2 + (^ (2) ) 2 + (3 



?(2)\2 



(A)n2 



A,a 



k(k-g)-A 
{k-g) 2 ■ 



When N — 3, they lead to a unique orthogonal solution 



£ 2 = (£ 



(2)\2 



(3)\2 



fc(A;-^)-4 



c(2) _ 



(3) 



(3) 



'2 J y - 3 ( fc _^)2 > -1 ""I " 

For iV > 3, it is convenient to introduce N — 2 dimensional vector notation 
/ E a 3) \ 

c(4) 

On 



0. 



(34) 
(35) 
(36) 
(37) 
(38) 

(39) 
(40) 



(N) 



(41) 
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The magnitude of £3 is given by 



l£sl ~ 3{k-g)> ■ (42) 



The magnitude of S 2 is controled by E\ through 

jk-g)- 
Kk-9) 2 



= k t: g) - A - (43) 



Similarly, we have, for 81, 

f+(g, w )'+igi'= y^ 4 . (44) 

They introduce three arbitrary constant parameters. We also know from 

£ 1 -S 3 = £ 2 -£ 3 = (45) 

that £3 is perpendicular to £\ and £2- The angle between £\ and £2 is fixed once we 
choose the three parameters controling their magnitudes by using 

£[ 2) £P + £ 1 -£ 2 = 0. (46) 

Apart from these geometrical conditions, we can take arbitrary combinations of the 
componenets for these three (N — 2 dimensional) vectors. In particular, there is 0(N — 2) 
invariance for the equilibrium values of £^\ Therefore, they are a family of equilibrium 
points that spans a 3(iV — 3) dimensional submanifold in the state space, for which 
the spatial slices are isotropic even though the vector fields have non-zero background 
values. 



3.2. Stability of the Isotropic Inflation 



To decide local stability of an equilibrium point requires linearized equations around it. 
Let us pick up one of the isotropic equilibrium points derived above. Before linearization, 
we shall manipulate the variables so that the stability analysis around that particular 
point simplifies. First, apply a rotation (-R) m ^ ne ^ ~ 2 dimensional internal space 



to introduce new variables T>^ 
satisfy 



(R)^^£a such that their values at that equilibrium 



( vf 





\ 



V 



J 



V 



(i) 



V 



J 








V 



(->) 







7^0. 



(47) 



/ 



This operation leaves the equations unchanged. When liearized around that equilibrium 



point, the perturbations for T>^ for B > 6 as well as T>^' and T>^' give trivially zero 
eigenvalues. Next, notice that at the equilibrium 

0. (48 



(3) 



p (2)p(2) 



p(3) p (3) 
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Thus we can choose a 3 by 3 orthogonal matrix Q such that 



1 c? 






( 


s[ 2) 


\ 




( 









) 




(s? 


\, 




( 





\ 






= Q 




P (3) 










MS) 




= q\ 










MS) 




U 4) 


J 




V 




J 




I 


) 




U 4) 


J 




v^ 4) 


J 




I 





/ 



(49) 



The linearized 



where the last equalities for each vector hold for the equilibrium, 
equations are invariant under this change and the eigenvalues for C[ A \ and 
are zero. Finally we take a 2 by 2 orthogonal matrix P to introduce a new pair of 
variables 

' £0-) 

£(2) 




(50) 



such that £^ = is satisfied at the equilibrium and its eigenvalue there is zero. 
After these transformations, we are left with £^\ Cf\ and T>a^ whose eigenvalues 
are potentially non-zero. All the other zero-eigenvalues are manifestation of the fact 
that we are considering a hypersurface of equilibrium points. Looking at the resulting 
linearized equations, we see that the problem is now reduced to the stability analysis of 
the orthogonal equilibrium point with three electric fields. The correspondance to the 
notation for iV = 3 is 

£V*>£, Cg^fg, Df'^f (51) 

The stability analysis for N = 3 is straightforward. The twelve independent 
varibales is eliminated by using the Hamiltonian constraint) are grouped into six 



J 12 



£1 



(2) 



J 13 



(3) 



pairs: 

- £ + £^ -i 
common eigenvalues given by 



J 2.3 



£ 



(3) 



£ 



(2) 



£ 



(3) 



2£ - £ 



(2) 



4 3) and $ 



£^ 3 \ The first five pairs, which represent anisotropic perturbations, share 



A 



S £ 



q - 2 ± y/(q - 2) 2 - 16£ 2 



Their real parts are negative since — 1 < q < 2. The isotropic mode has 

q - 2 ± y/(q - 2)2 - 4£^(4 - 3g(k - g)) 



A 



$ £ 



(52) 



(53) 



If we require q = < 0, it follows 4 — 3^(/c — g) > 0, which means that both of the 
eigenvalues have negative real part again. In other word, these isotropic equilibrium 
points are local sinks for a range of the parameters for which the universe undergoes 
rapidly accelerated expansion. 

One important consequence of this result (together with the global analysis in the 
next section) is that the anisotropy due to non-Abelian gauge fields should disappear 
if duration of inflation is sufficiently long since they can generically be regarded as a 
number of copies of Maxwell vector field, interacting with each other. The sufficient 
condition for isotropic final state is weakness of the Yang-Mills gauge couplings during 
inflation, which appears to be satisfied in the gauge- kinetic coupling model [22] . 
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4. The Cosmic Minimum-Hair Conjecture 

In this section, we would like to gather evidences for the cosmic minimum-hair conjecture 
which we would like to propose. In the first part, we clarify the global dynamical 
structure of the uniform coupling models. There, we see a mechanism of isotropization 
due to the generic structure of Einstein-Maxell equations. In the second part, we try to 
explore the fate of the anisotropy in models with more general couplings. We will see 
that the anisotropy of the universe is monotonically decreasing. 

4-1. Global phase space structure 

A natural question arises; whether those isotropic inflations with non-vanishing vector 
fields are the only attractor of the system. It is difficult to derive a global conclusion 
for a system of this complexity in general. It is also a parameter dependent problem. 
However, we can argue that multiple vector fields are expected to repel each other and 
try to become isotropic in expanding universes. 

Let us first look at the stability of axisymmetric inflating solutions which were 
discovered in [17]. In the present model, they are located on lower dimensional 
boundaries of the full state space. For example, we have an equilibirum point 

_5k 2 -2kg-3g 2 -8 2(k 2 - kg - 4) 12(fc - g) 

q ~ {k-Zg){k-g) + %' + ~ (k-Zg){k-g) + V ~ (k - 3g)(k - g) + 8' 

= 3((k + 3g)(k-g)-8)(g(k-g)-4) 3(k(k - g) - 4)((fc + 3g)(k - g) - 8) 

[(* - 3g){k -g) + 8] 2 ' [(k - 3g){k - g) + 8] 2 

£_ = E 12 = S 13 = S 23 = = S { a A) = 0. 

There are many others which represent physically the same spacetime, but lie on different 
boundaries. It was a stable attractor solution for the single- vector-field model. For this 
equilibrium state, we notice that S + > by looking at the evolution equation for S + 
and requiring S 2 > 0. More generally, any electric field in 1-direction tends to support 
positive S + by generating a tension along that direction. It is also easy to see that 
for any electric perturbation in 2- or 3-direction, the eigenvalue is given by 3S + > 0. 
Thus this axisymmetric solution is a saddle point in the general multi-field state space. 
Putting these together, the mathematical structure of the Maxwell's equations is such 
that the positive S + created by the 1-component of a field comes with negative sign in the 
eveolution equation in 1-direction while it has plus signs in the other directions. Thus, 
if there are more than one vector fields and one of them is dominant, it creates such an 
anistoropy that destabilizes orthogonal components of the other fields. This instability 
does not show up for single field models since it would merely cause a rotation of that 
vector. If those orthogonal components grow and surpass the original vector, it then 
becomes unstable and recovers its amplitude. 

Thus we expect that multiple vector fields always try to rearrange their orientations 
so as to minimize anisotropy of the space which would cause an instability in some 
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of their components. This is a generic feature of the Maxwell's equations in general 
relativity and has nothing to do with any couplings to other matter components. The 
point is that the anisotropy created by a tension acts rather to reduce the tension 
itself. In the present model, the vectors cannot disappear since they keep excited by 
the dominant scalar field. However, they still redistribute themselves to achieve as 
much isotropy as possible within the given circumstance. As an example, if N = 2, the 
attractor solution contains two orthogonal electric fields with the same amplitude: 

2k 2 - 2kg - 3g 2 - 2 k{k - g) - 4 3k{k - g) - 12 



<1 



(k -3g)(k-g) + 2' 2(k - 3g)(k - g) + 4' 2(k - 3g)(k - g) + 4 : 

$ = 6(fc - 2g) ^ _ 3(g(k - g) - 2)(g(2k - 3g) - 2) 



£ 2 = (£ 



(k - 3g)(k -g)-2' ((k - 3g)(k - g) + 2? 

{2) . 2 _ 3(k(k-g)-4)(g(2k-3g)-2) 



2((k-3g)(k-g) + 2y 



E12 — S13 — E23 — &i — 0. 

If we include a third vector field, however, it becomes unstable because of the eigenvalue 
for f 3 (3) 

A £s = 6E + > 0. (54) 

From the above argument, we can also read off the global behaviour of the orbits. 
There are a bunch of saddle points consisting of those axisymmetric signle-field and 
orthogonal two-field equilibrium solutions on the boundaries, which are attractors when 
restricted in properly chosen invariant subsystems. Since these invariant sets have 
dimensions smaller than 3(iV + 1), the initial condition from which the orbits are 
attracted to any of the anisotropic attractors is of measure zero in the entire state 
space. Whenever an orbit starts from a point not included in these subsystems, it is 
attracted to the isotropic final state after a sufficiently long time. For example, in order 
for an orbit to be attracted to the two field attractor, it must satisfy = for all 
A — 3, ■ ■ ■ , N. Any small deviation from this subsystem would lead it to isotropy due 
to the linear instability demonstrated above. A typical orbit is first attracted towards a 
nearby saddle point. Then the instability explained above kicks in and another vector 
component rises. It might come across another instability and go to another saddle 
point. An orbit continues this routine until it finally settles down to one of the isotropic 
attractors. 



4-2. The cosmic minimum-hair conjecture 

In this subsection, we numerically solve the dynamical equations fTT6l) - (1291) in the case 
of three vector fields. The results show that the cosmic minimum-hair conjecture indeed 
holds in many cases. 

First of all, we confirm the dynamical mechanism discussed in the previous section. 
We solved the basic equations with the parameters k = —2, g = 5, and the initial 
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conditions E± = E 12 = £13 = £ 23 = 0, $ = 0.6,5 = 0.1,£ x (2) = 0.2 ,£ 2 (2) = 0.001 ,e[ 3) = 

(2) (2) 

0.1 ,82 = 0.01 , £3 = 0.0001. The initial condition for \l/ is determined by the 
constraint equation. The parameters are chosen so that the universe is accelerating, 
namely q < 0. 




Figure 1. For a uniform coupling g = 5, we plotted time evolutions of anisotropic 
expansion rate normalized by the Hubble. The horizontal axis is the e-folding number. 
The components £+, £_, £23, £13, and £12 correspond to thick, dashed, blue, red, 
and green lines, respectively. 

From FigJTJ we see the anisotropy disappears after transient anisotropic inflationary 
phases. This confirms that the isotropic inflation is an attractor. The linear stability 
of the isotropic solutions shown in the previous subsection guarantees that the final 
isotropic plateau is not transient and never decays. It should be noted that the vector 
fields possess a nontrivial orthogonal configuration. We have solved the basic equations 
for other sets of parameters and initial conditions. Although the transient behavior 
depends on the parameters and the initial conditions, the system always achieves 
the isotropic inflationary final state asymptotically. The transient anisotropic phases 
correspond to the saddle points discussed in the previous section. However, since 
duration of anisotropic inflation is sufficiently long, the anisotropy at each of the saddle 
points would be relevant to CMB observations in the realistic cases where inflation ends 
with a finite duration. 

In order to see what happens for general couplings, we took couplings g\ = 4.8 , #2 = 
5.0 ,#3 = 5.2 for the three vector fields. As you can see in Figj2j the attractor is an 
anisotropic inflation. Strictly speaking, the final plateau in Fig J2] might decay eventually 
if we waited long enough since we didn't prove the existence of a stable fixed point for 
nonuniform couplings, although numerical evidence suggests otherwise. In this general 
case, the behavior is even more interesting. Since there are several saddle points, we have 
different anisotropic inflationary stages. During inflation the preferred direction in the 
universe is changing, which would lead to interesting observational consequences. For 
example, the anisotropy can be scale dependent. Even if we change initial conditions, the 
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Figure 2. The time evolution for the couplings gi = 4.8 ,172 = 5.0 ,93 = 5.2. The 
horizontal axis is the e-folding number. The components E+, S_, E23, E13, and E12 
correspond to thick, dashed, blue, red, and green lines, respectively. 



id 
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Figure 3. The time evolution of the total anisotropy E a& E a b/2 is depicted with respect 
to the e-folding number. 



attractor is the same. Hence, the predictive power of the theory is maintained. In Figj3j 
we plotted the square of the magnitude of anisotropy E af> E ab /2. While the system is not 
close to any of the fixed points, the orbit follows a complicated dynamical trajectory. 
However, once it starts to be attracted towards one of the inflating solutions, the 
evolution is fairly regular. The result shows monotonic decrease of the total anisotropy 
after the orbit goes into this regular tracking regime in spite that each component of 
the anisotropy exhibits a complicated behavior. It is inferred that when an orbit goes 
from a fixed point to another, the latter is always less anisotropic than the former. We 
have examined a lot of other parameters and initial conditions and confirmed that the 
total anisotropy is always monotonically decreasing at late times. 

The above evidence strongly supports the cosmic minimum-hair conjecture: that 
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is, the universe organizes itself so that any feature of the spacetime during inflation 
becomes minimum. The cosmic no-hair theorem can be regarded as a special case of 
this general principle. 

5. Conclusion 

We investigated inflationary models with multiple vector fields. In the presence of non- 
trivial couplings between the inflaton and the vector fields, the no-hair theorem does 
not hold any more. The universe can possess vectorial hairs in general. Thus, for a 
single or two vector fields, there remains anisotropy of the expansion. Remarkably, 
the expansion of the universe eventually becomes isotropic when the number of vector 
fields is greater than two and the couplings are uniform. More precisely, while one or 
two vector fields dominate the energy density of the vector sector, the expansion of 
the universe is anisotropic. However, since these phases are transient in the presence 
of other vectors, the trajectory eventually approaches the attractor, which has been 
shown to be the isotropic multi-vector inflation. This isotropization mechanism tells 
us that orthogonal triad is a stable configuration. This also tells us non-Abelian gauge 
fields do not create any anisotropy in the asymptotic future unless they remain in a 
strong-coupling regime during inflation. Curiously, a similar mechanism works in a 
different vector inflation model [32]. In view of these results, one may think that the 
no-hair conjecture survives. However, it should be stressed that there exists vectorial 
hair. The point is that vector fields isotropize the universe by themselves. In the 
case of many vector fields, the configuration cannot be an orthogonal system, rather 
the eventual configuration for vector fields depends on the initial conditions. In a 
sense, this is a "moduli space" for isotropic inflation with multi-vector fields. When 
the coupling constant g m for each vector field is different, we have anisotropic inflation 
as the attractor. Even in this case, there are several transient anisotropic inflationary 
phases. It seems that the inflating system dynamically selects a minimally anisotropic 
configuration. This has been confirmed numerically. On the ground of these analyses, 
we can expect that the cosmic minimum-hair conjecture holds. 

There is a caveat. In this paper, we have considered a power law type inflation, 
which is eternal. However, in the realistic inflationary models there exists the end 
of the inflation. Although the analytical treatment in the present article does not 
work for general scalar potential, the system would still be well-behaved and numerical 
calculation should be readily implemented. To study the end of inflation will be an 
interesting direction for further research not least because of its importance regarding 
reheating. We would observe anisotropy if the saddle solution earns a sufficient e- 
folding number and the inflation ends before the relevant scales go beyond the observable 
range. The phenomenological consequence of the anisotropy is worth being mentioned. 
It produces the statistical anisotropy and the cross correlation between curvature 
perturbations and primordial gravitational waves [251 123 ISU [28] • It would be interesting 
to see these features in the CMB [29, 33]. 
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